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Abstract 

We investigate the Hubbard model in the hmit U = oo, which is equivalent to the statistical 
condition of exclusion of double occupancy. We solve this problem using Dirac's method for con- 
straints. The constraints are solved within the Bosonization method. We find that the constraints 
modify the anomalous commutator. 

We apply this theory to quantum wires at finite temperatures where the Hubbard interaction is 
U = oo. We find that the anomalous commutator induced by the constraints gives rise to the 0.7 
anomalous conductance. 
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1. Introduction 



Advances in the physics of electronic devices requires the computation of the wave func- 
tion for confined geometries and strong interactions. A possible way to study such cases is 
the use of quantum constraints. A particular situation occurs in quantum wires when the 
short range electron-electron interaction is governed by a large repulsive Hubbard interac- 
tion. When this interaction obeys f/ — t- oo, double occupancy is prohibited. This problem, 



known in the literature as the Hubbard model at U 
method of Dirac's constraints 



oo 



2|-|8|, can be studied using the 



Interesti ngly the physics of quantum wire seems to depend on the electron-electron in- 
teractions {9131. The physics of electron-electron interactions in one- dimensional metals 
is described by a Luttinger liquid which replaces the traditional Fermi liquid description. 
For large Hubbard interactions, [/ — t- 00 at finite temperatures, one obtains an incoherent 



Luttinger liquid 



19 



20| characterized by magnetic excitations which are negligible in com- 



parison to the temperature. 

The physics of the Hubbard model at f/ — t- cxd might be relevant for studies of the 0.7 



conductance anomaly discovered in quantum wires 14l-ll7l|. Experiments suggest that the 



anomalous conductance is observed at low electronic densities 



18| and short in wires 24 1. 



Recent Monte Carlo simulations, performed in the limit of large Hubbard U and finite 



2lj. In ref.|22j. 



temperature, show that the one-dimensional conductance is anomalous 
the one-dimensional itinerant electron model with ferromagnetic coupling has also been 



investigated. Performing a Monte-Carlo study, the authors in ref. 22|] have found that the 
conductance in this case is anomalous as well . 

From the work of ref. 23|, we know that a one dimensional model with nearest neighbor 



hopping elements and [/ — > 00 has a ferromagnetic ground state at zero temperature. The 
strongly interacting Hubbard model U = 00 and the ferromagnetic interactions are both 
characterized by the exclusion of double occupancy. Therefore, it seems that the anomalous 
conductance reported by 21| and 22| might have the same origin. 



Experimentally, the anomalous conductance was observed in GaAs/AlGaAs at 



ow elec- 



25-27 



29|. 



tronic density where the long range Coulomb interaction plays a significant role 
In ref. 29(|, it has been shown that the long range Coulomb interactions in the limit of low 



electronic density drives the one dimensional wire to the [/ — ?■ 00 limit. 
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Inspired by the Monte-Carlo results 2l|, |22|], we will compute the conductance for a one 



dimensional wire in the limit f/ = oo at finite temperature using the method of Dirac's 
constraints. The limit U = oo gives rise to exclusion of double occupancy (this means that 
in the ground state the fermion occupation number takes only two values , zero or one ). 
Therefore, the statistics is similar to the statistics of spinless electrons. The exclusion of 
double occupancy is taken into consideration, by demanding that the ground state wave 
function \F > must be annihilated by the pair operator '$i{x)'^^{x)\F >= ( "^aix) is 
the single electron operator). Following Dirac 0, IsO], we learn that the constraints must 
be satisfied at any time. As a result, one generates a set of constraints which satisfy the 
equation Qi{x)\F >= ,i = 1,2,3.. 

In the language of quantum mechanics this means that one has to find the wave function 
\F > which is annihilated by the set of constraints Qi{x)\F >= 0, i = 1,2,3.. and is an 
eigenstate of the Hamiltonian H\F >= E\F >. When the commutator of the constraints 
[Qi{x),Qj{x')] with i ^ i can be inverted, the commutation relations are modified to the 



Dirac commutators 



301 ]. As a result, the Heisenberg equation of motion for an observable 



A is modified from = [A, H] to iJi^ = [A, Hjo where [A, Hj^ is the Dirac commutator 
given by [A,H]d = [A, H] - [A,Q,]{[Q,,Q2])-'[Q2, H] - [A,Q2]{[Q2,Qi])-'[Qi, H]. The 
presence of the constraints [Qi, Q2\ (in the Dirac commutator) modifies the Pauli statistics 
to the spinless electrons statistics. 

For U = oo, the Luttinger liquid at finite temperature is replaced by the incoherent Lut- 
tinger liquid characterized by the pinned spin excitations and propagating spinless electron. 
Therefore, only the electron current is conserved. As a result, when an electron with spin 
up or spin down is injected at one end of the wire, we find that the electron which reaches 
the other end of the wire has an arbitrary spin. The interactions are restricted to the wire 
and are absent in the leads. Therefore, the density in the wire corresponds to the spinless 
electron density which is half the density of the non-interacting electrons in the leads. This 
change of the density will modify the commutation equation for the Bosonic densities by 
a factor of two. The effect of this modification is observed at a finite temperature when 
a voltage difference is applied between the two reservoirs which are connected to the wire- 
leads system. The injected current from the reservoirs will be twice the transmitted current 
into the wire. The interface between the leads and wire randomizes the spin and causes the 
electrons to behave as spinless particles at finite temperatures in the wire. Therefore, the 
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conductance will be half in comparison with a non -interacting wire. 

The plan of this paper is as follows: For pedagogical reasons, we present in chapter 2 the 
method of Dirac's constraints used in quantum mechanics. This method will be used in the 
remaining chapters. In chapter 3, we present the model for an interacting wire {U = oo) 
of length d coupled to the non- interacting leads and reservoir [U = 0). In chapter 4, we 
derive the anomalous Dirac commutator for the quantum wire model. This is done by 
extending the Bosonization method to the constraints . In chapter 5, we identify the electric 
current operator and in chapter 6, we compute the static current at a finite temperature. 
In chapter 7, we consider the effect of the Zeeman interaction, and in chapter 8, we present 
the conclusion. 

2. Dirac's method in quantum mechanics 

In this section, we will present Dirac's method in quantum mechanics. This method will 
be extended in the next sections to the strongly interacting quantum wires. 

We consider the Hamiltonian H with the constraint operator Qi. The constraint op- 
erator is defined by the condition that when it acts on the true ground state |F > it 
satisfies the equation Qi\F >— 0. Formally this condition is introduced by demanding 
that the ground state wave function |F > of the Hamiltonian H has to obey, in addi- 
tion the constraint equation, Qi\F >= . Since the constraint must be obeyed at any 
time, we must have ^!:9lMl^p >= q. From the Hcisenbcrg equation of motion, wc ob- 
tain the condition ih^^^^^^\F >= [Qi{t), H]\F >= 0. Therefore, wc find that the con- 
straint is satisfied at any time, only if the constraint operator Qi commutes with the 
Hamiltonian [Qi{t),H] or satisfies the equation [Qi{t),H] oc Qi{t). If this is the case, 
Qi will be the only constraint. Here we will assume that this not the case! There- 
fore, the only way to satisfy the condition ^^i^-^F >= o is to introduce a second con- 
straint Q2 oc [Qi,H] (or Q2 oc [Qi,H] — constant ■ Qi )■ For simplicity we assume that 
no more constraints are needed. As a result, we obtain a problem with two constraints 
Qi\F >— 0, and Q2\F >— which do not commute [Qi,Q2] 7^ 0. To enforce the 
constraints on the wave function, we use the method of Lagrange multipliers Ai and A2. 
This allows us to replace the Hamiltonian H by Ht — H + XiQi + \2Q2- 
We will determine the Lagrange multipliers using the Heisenberg equations of motion. 
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Prom ih^\F >= ([Qi,//] + A2[Qi,Q2])|i^ >= 0, we obtain A2 = -{[Qi,Q2])-^[Qi, H], 
andbomih^lF >= {[Q2, H]+Xi[Qi,Q2])\F >= 0, we obtain A2 = -{[Qi,Q2\)-^[Qi, H]. 
We substitute the explicit form of the Lagrange multipUers into the Hamiltonian Ht, and 
find the following Heisenberg equation of motion for any operator A: 
dA 

ih— = [A,Ht] = [A,H]-[A,Q^]{[Q,,Q2])-'[Q2,H]-[A,Q2]{[Q2,Qi]r'[Qi,H] = [A,H]n 

(1) 

The conclusion from this calculation is that the effect of the constraints have changed the 
regular commutator [A, i^] to the Dirac bracket commutator [A,H]iy, which is given by 
[A,H]n = [A,H] - [A,Q^]{[Q,,Q2])-'[Q2,H] - [A,Q2]{[Q2,Qi])-'[QuH]. 

Using this result, we define the Dirac bracket for two operators A and B in the following 
way: 

[A,B]d = [AB] - [A,Q,]{[QuQ2]r'[Q2,B] - [A,Q2]{[Q2,Qi])-'[QuB] (2) 

This methodology will be applied in the next sections to the problem of strongly inter- 
acting electrons. 

3. The model for a short wire with exclusion of double occupancy coupled to 
non-interacting leads and reservoir 

We consider a quantum wire of length d which is perfectly coupled to the leads of length 
L — d shown in figure 1. The leads are in thermal equilibrium with the two electronic 
reservoirs which have chemical potentials fiR and fi^. The reservoirs arc described by three 
dimensional non-interacting electrons. Physically, once the electrons reach the reservoir they 
do not return to the wire . At a finite temperature T, the infinite leads can be replaced by 
finite leads L ^ Lt {Lt ''iS the thermal length for which coherency in the leads is preserved). 
At finite temperatures we have the model: H — Hy^ire + Hig^ds + -f^r 



' reservoir 



(3) 

U{x) is the space dependent Hubbard interaction, n^ix) — '$]^{x)'i>aix) represents 
the electronic density for the spin polarization a —^,-1 and Ep is the Fermi energy. 



The Hubbard interaction U{x) is restricted to the region |x| < | (with the condition 
d « L). In the comphmentary region |x| > | the Hubbard interaction is U{x) = 0. 

In the remaining part of this section we will present the Bosonized version of two models 
which have the same Hamiltonian and differ only by the constraint conditions. The two 
models are: the non interacting electrons described by [/ = and the strongly interacting 
electrons U = oo. Both problems will be investigated using the method of Bosonization. 
The method of Bosonization is based on the representation of the Fermion operator \E'cr(a:) 
in terms of the Bosonic fields '&Rn_= ^'^7'^'^ , 'dia = ^'^t'^'^ ■ The Bosonic representation of 



the electron operator is given by [3ll433l| : 



where a is the lattice constant and Kp is the Fermi momentum. Using the Bosonic fields 
'^L,a{x) and t?i?,o-(a^), we define a new field Pa{x) = dx(pa{x). In order to construct a quantum 
theory, one has to identify the canonical momentum operator which is conjugate to the 
Bosonic fields '!?o-(x). Here, we have an example where both the operator ^cr{x) and Pcr{x) = 
dxfa{x) (the candidate for conjugate momentum) are built from the same fields ( t^ij^o- ,'*?l,(t 
). Therefore, the commutator [d^lx), P„/{x')] can be finite only for an infinite number of 



particles and vanishes otherwise. Such commutators are called anomalous [31|, 133|]. For such 
cases the commutator ['i9o-(x), Pa'{x')] will be defined according to the many particle ground 
state \G > with the expectation value < Glldaix), Pa'{x')]\G >. 

A-The non-interacting electrons, U{x) = 

The Hamiltonian for the non-interacting electrons is given by eq.3 with U{x) = 0. The 
Bosonized Hamiltonian is given by : 

where tlrA^^ ^^'^ ni^^^x) are the right and left electronic densities and ///j, yWL are the 
chemical potentials. The non-interacting electrons ground state is identified with the Fermi 
sea \F^^^ >. We compute the ex pec tation value of the operators [d^{x), Pa-'{x')] and identify 



the anomalous commutator 
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33|: 



< fW|K(x),P.,(x')]|FW >=i6^yS{x-x') 
6 



(6) 



The result in equation (6) allows to identify the operator Pcr{x) = dxfcrix) as the canonical 
conjugated momenta operator and define the commutator [d^^x), Pa'{x')] = i6a,a'S{x — x'). 
(Once the commutators have been defined we can choose the representation representation: 
'd^{x)\'d^,'i!}^ >= i3^(x)\'d^,'di > with the conjugated momenta given by Po-{x) = —ij^^j^-) 

B-The strongly interacting electrons U{x) = oo 

We consider the strongly interacting case described by the Hubbard interaction U (x) = 
oo. ( For finite Hubbard U interactions, the Coulomb interactions drive the the Hubbard 
U interaction to U — > oo. The limit U{x) = oo describes effectively the physics in the strong 
coupling limit ) 

The Hubbard interaction U{x) = oo is restricted to the region < |. As a result the 
the electron occupation number is restricted to the charge values ne{x) = 0,1. 

^U(^)^-t(^) + '^i=i{x)^a=dx)\F >= n,{x)\F > (7) 

Where |F > is the electronic ground state which is different from the non- interacting ground 
state \F^^'^ >. This condition is named the exclusion of double occupancy and is implemented 
with the help of the constraint equation. 

d 

'^l{x)'$^{x)\F >= for the region \x\ < - (8) 

In the complimentary region > | the electrons are non-interacting, and we have 
U{x) = 0. 

Once the constraint conditions have been introduced we can use the non-interacting 
Hamiltonian for the wire-leads system. 

H^^re + Hi^^ds ~ l\ dx[ ^ ll^ml{x)^^> „{x) - Ep (9) 

Next, we Bosonize this model using the same fields as we used for the noninteracting 
case. In spite of the fact that the Hamiltonian for the non-interacting case is identical to the 
strongly interacting one, the two ground states are different. The new ground state |F > 
obeys the constraint equation (8) which, in the Bosonic form, is given by the representation 
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^^{x)'^l{x)\F >= 2e-'^''^^''^cos[2KFX + v^t9e(x)] + cos[V2^'d s{x)]\F >= 0; \x\ < ^ 

(10) 

The constraint operator is represented in terms of the Bosonic fields for the charge degrees 
of freedom: i?e(x) = {^^{x) + ^i{x)) , Lpe{x) = l{y^t{x) + ipi{x)) and Pe{x) = d^i^eix). 
Similarly, the spin degrees of freedom are given by: ^s{x) = ■^{'^^{x) — 'di{x)), ips{x) = 
■^{ip^{x) - ipiix)) and Ps{x) = d^ips{x). 

In order to find the commutation relation, we will use the new ground state which is 
restricted by the constraint. The commutator [drrix), Prji(x ' )] w ill be defined by the value 



of the expectation value with respect the ground state [3ll, |33|]. This expectation value is 
the central result of this paper and is derived in the next section. The expectation value 
of the commutator is a function of h{x,d) given by: h{x,d) = 1 for |x| < ^, {U = oo) and 
h{x,d) = for |a;| > |, (?7 = 0). 

< F| [Mx), Pe{x')]\F >= i6{x - x')[l - h{x, d)^] (11) 

This result allows us to define the Dirac commutator [deix), Pe{x')]Dirac = iSi^x — x')[l — 
h.{x, d)^]. 

The Bosonized Hamiltonian for this case is given by : 

H^^re= E f] dx^[{d,M^))' + {d,M^)f] (12) 

(13) 

V is the Fermi velocity for the non-interacting fermions. We assume perfect coupling between 
the wire and the leads. The Bosonic fields obey continuous boundary conditions: da^x = 
f-e) = d^{x = f + e), 'd^{x = -f-e) = ^a{x = -f + e) and y?<,(x = f-e) = ip^{x = f + e), 
(Pa-{x = — e) = (faix = + e) , where e represents the overlapping wire- leads region. 
Therefore, we can replace the Hamiltonian Hy^i^e + Hi^ads by a free Bosonic model. 

^ hv , 



H^ire + Hieads = E /! dx-^[{d^Lp„{x)f + {d^,d„{x)f] 

Hreservoir= E j _^ dx[HL{,x)nL,a{,x) + ^R{x)nR^„{x)] (14) 



Hreservoir IS the reseivoir Hamiltonian which is a function of the voltage difference 
ij'r{x)-^.l{x) _ y right and left electronic densities nji^^{x) = PR^a{x)+ < 

F\nR^a{x)\F >, riL^f^i^x) = pL,a{x)+ < F\nL,a{x)\F >. The Bosonic representation of the 
densities is: pR,aix) = j^{d,^'d^{x) - d^ip^{x)) and = ^{d^-d^ix) + d^ip„{x)) 
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33| . Due to the fact that the commutator for the electronic densities is space depen- 



dent, the average density is space dependent < F\nRfj{x) + nL^fj{x)\F >= ne[l — ^^"^'^^ 
where ne = — is the electronic density in the leads. The space dependent electrostatic 



potential caused by the space dependent charge density is given by bv{x) = |3|. 
This potential will introduce current fluctuations. 

4. Dirac's method for the exclusion of double occupancy-an application to 
the finite wire-leads system-A Bosonization formulation for the constraints 

We will extend the Quantum Mechanical results obtained in the previous section to 
the wire Hamiltonian H^ire and we will introduce a new formulation for the constraints 
using the method of Bosonization. For the wire Hamiltonian Hy^ire, we will enforce the 
exclusion of double occupancy. The operator for exclusion of double occupancy js] leads to 
the constraint condition for the ground state |F > which replaces the non- interacting Fermi 
sea \F^^'^ >. The exclusion of double occupancy is enforced by the two electrons operator 
'^^{x)'^l{x)\F >= for l^l < I . Except for half filling, \E'f (x)\E'4,(x) is the only primary 
constraint Jj. In order to find all the secondary constraints [1], we have to commute the 
constraint operator '$^{x)'$i{x) with the Hamiltonian. If this commutator is not equal to 
the constraint field, new constraints are generated. Using the Bosonic representation \E'o-(a;) 
given in equation (5), we compute the representation of the pair operator '$^{x)'$i{x). The 
pair operator is represented in terms of the Bosonic fields, 'de{x) , ipe{x),'ds{x) and fs{x): 

'^^{x)'^lix)\F >= 2e-^^^'=(^)Qi(2:)|F >= 0; |x| < ^ (15) 

where Qi{x) is the Bosonic constraint : 

Qi{x) = cos[2Kfx + ^/^^e{x)] + cos[V2^'ds{x)]; Qi{x)\F >= (16) 

In order to understand the effect of the constraints, we will present a simplified description. 
For this purpose we will ignore the 2Kf oscillations, and we will approximate Qi{x) by 
Qi{x) ~ cos[\/27r'ds{x)]. From the equation cos[\/27r'ds{x)\F >^ 0, we learn that the 



spinon phase must be pinned to a constant value '!9s(x) = (v27r)~^[| + nir], n = 

0,1,2... Therefore, the system has no spinon excitations, and the only excitations are spinless 
fermions. It is this constraint which gives strong exponential decays of the spinon correlation 
function. As a result, the system is described effectively by a spinless system. Therefore, 
due to the exclusion of double occupancy, the electronic density in the wire will be half 
in comparison to the leads, where the constraint is absent. This changes will modify the 
commutation equation for the electron number from [d^i^), Pe{x')] = i6{x — x') to a new 
commutator [&eix), Pe{x')]D = — x'). 

For the remaining part, we will work with the full constraint operator given in equa- 
tion (16). The constraint Qi{x,t) must be satisfied at any time. Therefore, we must have 
>= 0. Using the Heisenberg equations of motion with the non-interacting anoma- 
lous commutator given in equation (6), we compute : ih ^^^j^'^^ = [H^ire,Qi{x,t)]. The 
condition ^Q2lEA\F >= can be satisfied only if one introduces a new constraint Q2{x) 
given by [Qi{x,t), H^ire] ■ 

Q^(x) = V2sin[2KFX + ^/^^e{x)]Pe{x) + sin[V2^^six)]Psix) (17) 



The constraint must annihilate the ground state, Q2{x)\F >= 0, and we must also have 
>= 0. In order for this to happen, we need to include a third constraint Q3{x) oc 
[[Qi{x,t) , H^ire], H^ire] which is given by: Qsix) = sin[2KFX + y/T:de{x)]{y/2Pe{x)Y + 
sin[\/2Ti'ds{x)]P'^{x). Continuing this process, we identify additional constraints Qa{x)^ 

Q^i^x) The third order constraint Qz{x,t) oc [[Qi{x,t), Hyji^f], Hwire\ can be written as 

a second order time derivative Qz{x,t) oc LQi^^_ xhe higher order constraints Qi{x,t), 
i = 3,4, 5, .. can be neglected since they are represented by higher order time derivatives : 

Following [l|, |30|, we find that the operators Qi{x) i=l,2 form a set of Second Class 
constraints, [Qi{x),Q2{x)] ~< F\[Qi{x),Q2{x)]\F >^ 0. In order to find the new commu- 
tation rules for our system, we will use the Lagrange multiplier fields Aj(x) which will be 
multiplied by the step function \\{x, d) (which is one for |x| < | and zero otherwise). We are 
interested to find the commutation rules for the full wire-leads system. In order to do this, 
we replace the Hamiltonian H^ire + Hi^ads by the the Hamiltonian Hx which is a function of 
the Lagrange fields Aj(x) that enforce the constraints Qi{x) for the wire |l|. Since we have 
no constraints for the leads, we will introduce the step function h{x, d) which is zero for the 
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leads and takes the value of one for the wire. The constraints for the wire will be enforced 
by replacing the Lagrange fields \i{x) by the product h(a:, d)\i{x). 



, L 



Ht = H^ire + Hi^ads + / '^ dxh{x, d) ^ \{x)Qi{x)) (18) 

''~ i=l 

We will apply the quantum mechanical method given in equation (9) to the Hamiltonian 
for the wire given in equation (18). Since the constraint must obey f^^i^hlij^ g, and the 
commutator [Qi[x),Q2{x')] can be inverted, the Lagrange multipliers can be computed. 

.^ dQi{x,t) ^^ >= [Q,{x,t),H^,re + Hieads + dxh{x , d)^. K{x)Q,{xmF >= 

(19) 

Prom equation (19), we obtain the Lagrange fields Xi{x). The Lagrange fields are functions 
of the matrix Ci^r{x, u), r = 1,2 which are defined with the help of the the integral equation 
J2l=i J duCi^ri^, u)[Qr{u), Qj{y)] = 6ij6{x—y) ,j — 1, 2. The matrix Ci^r{x, u) is given by the 
inverse of the commutator [Qi{x), Qr{x')\. The commutator of the constraints [Qi{x), Q2{x')\ 
(given in equations (16 — 17)) is computed using the non-interacting Bosonic commutator 

[Qi{x),Q2{x')] = -i^/2^8{x - x')[{sin{^p2^^s{s))f + {sin{2KFX + 0F^?e(^)))'] (20) 

We substitute in equation (19) the result of the commutators [Qi{x), Q2{x')\ given in equa- 
tion (20). As a result, the equation i = 1,2 ( see equation (19) ) allows us to 
determine the Lagrange fields Xi{x), i = 1,2. Next, we substitute the Lagrange fields Xi{x) 
into the total Hamiltonian Ht given in equation (18) and find that the Heisenberg equa- 
tions of motion are determined by the physical Hamiltonian Hyjire + Hi^ads with the modified 
commutation relation [, Jciroc- The Heisenberg equation of motion for any operator A{x,t) 
will be given by : 

[A{x,t),H^ire + Hieads]- J2 J2 / duh{u,d) dv[A{x ,t) , Qi{u)]Cij{u, v)[Qj{v) , Hy,ire + Huads] 

i=l,2j=l,2'' •' 

(21) 

The Dirac commutator for the canonical conjugates fields 'i?o-(a;) and Pcr'{y) wiU be given by: 

[M^),Pa'{v)]Dirac = K (x) , P,. (y)] - f dvh{u,d) [_[ dv[Mx) , Q2{u)]C2,l{u, v)[Qi{v) , PAv)] 

(22) 
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We observe that the Dirac commutator given in equation (21) has exactly the same structure 
as equation (2) obtained in chapter 2. 

The exphcit form of the Dirac commutator for the electronic density operator 'de{x) and 
the conjugate field P(.{x'), [deix), Pe{x')]Dirac will be given by: 

h{x,d){sin{2KFX + ^/7^^e{x)))'^ 



(sm(v^T9,(x)))2 + {sin{2KFX + ^de{x))Y 

(23) 



The constraint operator Qi{x)\F >= causes any operator R{x) to obey the equation: 
R{x)Qi{x)\F >= 0. We choose R{x) = cos[2Kpx + ^/7rl9e{x)] — cos[\/27ri!)s{x)] and con- 
struct the product R{x)Qi{x) = {cos[2Kpx + y/7T{}e{x)])'^ — {cos[y/27r'ds{x)]y . Adding and 
subtracting one from R{x)Qi{x) gives us a new equation R2{x) which obeys: 



R2{x)\F >= [{siniV27T^s{x))f - {sin{2KFX + ^/T^{}e{x))f]\F >= (24) 

This gives as the result: 

{sin{2KFX + V^i^e{x))y ^.^ 



{sin{V2^{}s{xW + {sin{2KFX + y/^de{x))f' 

{sin{2KFX + ^/^^e{x))y , ^ 1 

^ -L ^ 

Of QonOK^T -I- .^')') ('r^^^'i\^ l (isin{V^^s{xW-isin{2Kpx+y/^^e{xW y 2 



(25) 



The result in equation (25) was obtained after expanding the last expression in pow- 
ers of ( (sin(2Ki!x+^§ (x)))^ ] )" using the condition {R2{x))^\F >= The commutator 
[deix), Pe{x')] is finite only for an infinite number of particles and vanishes otherwise. The 
anomalous commutator will be given by the expectation value < FU'deix), Pe{x')]\F > 
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33| . This allows us to define the Dirac commutator [deix), Pe{x')]£)i 



[M^),Pe{x')]Drrac =< F\[d,ix) , P,ix')]D^rac\F > 

/NT , / ,N ^1 (sin(2KFX + ^yni!}Jx)))'^ , 

= id(x - x') l - h(x, d) < F =^ ^ ^ ) F > 

^ ^ ' \sin{V2^^s{xW + {sin{2KFX + ^MxW'^ ^ 

= iS{x-x')[l~^h{x,d)] (26) 

This result shows that the anomalous commutator (Dirac commutator) for the wire has been 
modified to ^h. This is consistent with fact that the electronic density for the wire is half 
the density of the leads (see figure 1). 
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The Dirac commutator affect in a significant way the equation of motion for the particle 
-hole excitations. Using the Heisenberg equations of motion ih^^^^^^ = [di,{x,t), H^ire]Dirac 
and ih^^^^^^ = [(peix,t), H^ire]Dirac wc find: 



[d^-v\l-^^rdl]M^,t) = (27) 

The space dependent commutation rules give rise to a multi particle-hole scatter- 
ing state. Performing a space average, we obtain the effective velocity Veff. = 
f y< (1 - > spacer average = v{l - ^) with the fundamental frequency u = Veff.^,. 

The space dependent part in equation (27) gives rise to the scattering potential Vsc{x) = 
v'^[< (1 — >^p^^g_^^g^^^g —(1 — h(M) -|^2j ^j-^jf,];^ generates high harmonic states. 

To conclude the anomalous commutator for a wire with exclusion of double occupancy 
given in equation (26) is the central result of this paper. This result was obtained restricting 
the infinite number of generated constraints to a finite set of constraints Qi{x) and Q2{x). We 



have shown that according to scaling theory 35|], the higher order constraints are irrelevant 
and therefore can be neglected. 

5. The current operator for the wire 

Using the Dirac commutator, we find the equation of motion for the electronic density 

Pe{x) = -^d^^eix). 
dpe{x,t) 1 



dt ih 



a,( dylMx), Pe{y,t)]Dn'acPe{y,t)) = -d^Je{x,t) (28) 



We identify from the continuity equation ^^"j^'*-* + dxJeix,t) = the conserved electronic 
current operator Je(x). We multiply Je{x) by the electric charge (— e) and find the electric 
current operator /g. 

/e = ^ /* dy[Mx),Pe{y)UracPe{y) = ^[1 " ^^^]d,Mx) (29) 

This equation shows that the current in the wire is density dependent. We find that the 
current in the wire (|2;| < |) is given by le = ^^dxf^eix) is reduced by a factor of two in 
comparison with the non-interacting electrons (|x| > |) 1^'^=°) = -^dxi^eix). 
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Contrary to the non-interacting wire where both the charge and spin currents are con- 
served, the constraints in the wire allow only for the conservation of the electric current Jg. 
Due to the fact that the commutation equations and the conservation laws for the currents 
in the leads and in the wire are not the same, we have at the interface between the wire and 
the leads incoherent scattering. (Electrons with a well defined spin (spin up or spin down), 
are injected from the reservoirs. The leads-wire interface randomizes the spin and causes 
the electrons to behave as spinless particles at finite temperatures.) 

6. The current at finite temperatures for a finite wire L » d and L > Lt 

In this section we consider explicitly the limit of f/ — oo taken at finite temperature 
T. In this limit, spin exchange processes are suppressed and the only allowed processes 
are the particles-holes (non-zero modes ) and particles excitations (zero modes). Due to 
the fact that the system is finite, we can separate the Hamiltonian, and the Bosonic fields 
into the periodic non-zero modes Bosons ^r^'^\x) = '&R,a{x), = '^L,cr{x) and the 

zero mode Bosons 'd^^~'^\x) = ^r^A^ij.cr, (x) = ^^^N^^^ where Nji^^j N^^^ measure the 



added charges with respect the Fermi energy [32|, l33 |. 



Using the Boson decomposition we find the representation: 



d,^e{x) = a,v^("^°)(^) + d,^lr='>\x) = d,0e{x) + ^[Nt,L + N^,L - N^,R - N^,r] (30) 
Similarly we decompose the Hamiltonian into two parts: H^j^^'^ + h\^^^^ = 

2 

and the zero mode part H^^i;^^ + H^'^^J^J: 



Hl^::^ + <r = ^ E + Nln] + ^^K, + Nl^]) (31) 

o-=t,4- 

For a constant external voltage, V = we have only the zero mode reservoir 

Hamiltonian: 

HiT^lir = E [f'LN^^L + ^rK,r] = E [i^)N.,R - i^)N^,L] (32) 

Since the commutator is space dependent, we will introduce a space dependent electrostatic 
potential 5v{x) = ^^^^^H-/^-^^^) which will couple to the non-zero mode densities. The effect 
of the electrostatic potential 6v{x) is introduced using the non-zero mode reservoir : 

Hiltervair = j\ dx[{-e)6v{x, t){ne{x, t) - 71^)] (33) 
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The electronic density is defined with respect the density in the leads: ne{x,t) — Ue = 
-^dx'&eix,t) + ne[l — ^^"^'^-^ ] — Tie = -^d^de jx, t) — TT-e ^^"^'^-^ ■ The electrostatic potential 
6v{x,t) is determined by the Poisson equation [34| : 



2 



K is the dielectric constant and ^ is the gate screening length [27l |. 

Next we will compute the two components currents dx(p^J^^^\x) and d^^'^^'^^x). The 
non-zero mode part dx^p'^J^^'^^x) is induced by the non-uniforni ground state —n j^^^„^^^ . 



This current will be obtained within the linear response theory 
the regular commutator by the Dirac commutator : 



, where we have replaced 



6 « F\d:,0e{x,t))\F »T=« F\dx0e{x , t)\F »TM^,t) " « F\dx0e{x , t)\F »T 

i ri; 1 hfv d") 



2 

, L 



r'^ll - r dy^^dxl , ' - , ' ] (35) 



2 



<< >>T stands for thermodynamic expectation value. Using the low energy Bosonic 

spectrum, which emerges from equation (27), we approximate for times t obeying Veff.^t < 
1, the time integration by a constant time r if Vefj.^r < 1. At finite temperature, r is 
given by r = sma//est[— -^1 = — ^ The zero mode current is computed using the 
current operator dx^^^^'^\x). At finite temperatures static electric current will be given by 
the thermodynamic expectation value of the operator ^t.-c-+^-i-,£^^t.-R ^-i-r ^ This expectation 
value is determined by the thermal reservoirs. The thermal function is represented by 
<< F\{N(j,L — Ncr^R)\F »T (T stands for the thermal occupation values). For the wire 
region, we have the Hamiltonian -ff£"7e°^ = |f I]o-=t,4,(i)[^CT.L + ^an] "which is subjected to 
the constraint of exclusion of double occupancy. For the leads, we have the unconstrained 
Hamiltonian H^Zt^ = ^ j:.=u{^)[NIl + N'^,r]- In the limit of long leads | « 1, the 
thermal weight of the region d is negligible in comparison to the leads region. Therefore, 
the thermal expectation value is given by the unconstrained Hamiltonian: 

^(n=0) ^ ^ ^ [^2^ + ^2 ^] + J2 [f^LN.,L + /ii?iV.,ii] (36) 

We compute the partition function under the assumption that 4 << 1. Therefore, the 



^(n=0) (n=0) n = 0) 

leads reservoir 



constraint partition function Z = Tj.[e ] can be replaced by the un- 
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constrained partition function Z = Tr[e ^'s^ ] computed with the unconstrained Hamil- 
tonian given in equation (22). As a result, the thermal occupation function << 

F\{N^,L - N^,r)\F »t^ ^TrliK^L - N^,R)e 



H-("=0) 



is given in terms of the thermal 



Fermi-Dirac occupation functions fp.D.i^n — ^I'l) 



1 



fp.D.i^n — IJ-r) 



1+e 



expressed in terms of the single particle spectrum en 

N.,R 



^(n-i),n=l,2,3... 



32 



33|: 



I h("=0) 

F »T^ -^TriiK^L - K,R)e 



"=°° eV 2L 

^1 I^T h V 



(37) 



In order to compute the static current caused by the static potential V, we have to perform 
a space average which singles out the zero mode current. This can be achieved if we perform 
a space average over the length Lt,l- After performing the space average, we see that on 
a length scale Lt that particle-hole current I^'^^'^^ vanishes. This result is seen from the 
the result v}^^^ + vj:Jj^ = given in equation (38). Therefore, we find ( see equation (38)) 
that the thermal expectation value of the current is given by zero mode current operator 

Using the length Lt,l = smallest[L, Lt], we perform a space and thermodynamic average 

/ 



-^T,L Lt,L 
2 



« F|/e|F >>T and find: 
— dx 



dx 



« F\IJF »T 



I" -^«F|/("^°)|F»T+ -^«F|/("=°)|F»T 
7^ Ltj, ' ' ' Ltj, ' ' ' 



27rp 

(--)(^)(— ) 



V 



h 



e^Up 



K 

,27re. 



+(r— ) 



V. 



eff. 



h 



dx{l — 



dx{l 



h(x,d) 



h(x,d) 



h(y,d) 



2 ^ J{x — yY + a? 



h(y,d) 



h(y,d) 



h(y,d) 



^{x~yf + a^ ^{x-yf + e 



+- 



ev 



dx_^^_h{x^^ E «F\N.,L 



N^,r\F »t 



L 



dx h(x^^^2^y) = [1 



T,L 



2 h^' 



d 



(3J 
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V^jj and V^jj represent the opposite voltage drop on each side of the wire which, due to 
the symmetry h(x,d) = h(-x,d) obeys, Vg^j + V^jj = 0. 



^ne(^)(-)iy[§,^] (39) 

In the present case, due to the symmetry, we can not observe the voltage ^e(^7^)(^)W^[f, „]• 
We believe that for an asymmetric situation this voltage might be measured. For the case 
where half of the wire-leads system a; > | is controlled by the strong interaction U = 00 
and the other half x < | is non-interacting U = 0, the conductance will be affected by the 
voltage . This voltage is in the range of lO^^volts and might be measurable. In figure 
3 we show the potential ^] as a function of the screening length. 

Equation (38) shows that the conductance for — = 1 is given by 6" = 2(0.5^), reflecting 
the fact that the anomalous commutator for the wire has been modified to ^h. For different 



lengths the conductance is modified. Following |2l|, we investigate the dependence of the 
conductance in terms of ^ and (the effect of the overlapping region e between the leads 
and the wire studied by 21 1 will be neglected here). The thermal length Lt = -jf^ = 
= (1 - where lJ'^ is the thermal length defined in the leads. The 

length L-Ti ^ Lt < L is given by the crossover function: ^ -^[1 + As a result we 
find for the conductance: 



(40) 



The result, of equation (40) is shown in figure 2. We observe that when the length of the 
wire d = -j^^ is in the range ^ ~ 1 — 2 the conductance is anomalous in agreement with the 



experimental results reported in ref. [2J]. (In the experiments the anomalous conductance 
is obtained as a result of the variation of the gate voltage, we will interpret this result as the 
region where the strong coupling point U = 00 has been reached. At this point, when the 
temperature or the length are varied the conductance varies. We believe that this variation 
of the conductance is demonstrated in figure 2.) 
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7. The effect of the Zeeman magnetic field 



A complete discussion of the anomalous conductance must include the ferromagnetic 

n n 

interaction proposed by [22|. According to [22], the emergence of the 0.7 conductance 
anomaly is due to ferromagnetic interactions. Since the effect of constraints takes into 
consideration the exclusion of double occupancy, which also holds for the ferromagnetic 
case, we believe that our results include, in a qualitative way, the results obtained in ref. 



22| 



The results reported by [22| show that in the presence of Zeeman magnetic field the 0.7 
plateau evolves into a robust 0.5 plateau. In our case, the value of the conductance varies 
from 0.5 to 0.7 (in the absence of the Zeeman interaction). These values are determined by 
the function h(a;, d) the temperature and length of the wire. 

In order to consider the effect of a Zeeman interaction, we observe that the Fermi velocity 
V is replaced hy: = v + ^, = v - ^ where A = {^)Kf[\[iV'-^^^ - . /i _ 



represents the Zeeman interaction, Kp = ^f^^f is the Fermi momentum, B\\ is the magnetic 
field and irip is the electronic mass. As a result, the Hamiltonian for the wire takes the form: 



(41) 

This Hamiltonian must be accompanied by the constraint operator of exclusion of double 
occupancy. 

^-,.(2;)*4(a;)|F >= 2e-'^^^(^)gf ~(x)|F >= 0; ^<x<^ (42) 
where (5f^^™''^"(x) is the Bosonic constraint in the presence of the Zeeman interaction: 

■\Zeemanf\ _ „^„ror^ _ i . /o — Q f i ^^„r/ ^\A™ i . /o — Q /_M. r^Zeeman/ 



Qf~{x) = cos[2Kfx + V27r^,{x)] + cos[{^)Ax + V2n^,{x)]; Qf~{x)\F >= 

(43) 

The Hamiltonian and the constraint operator show that neither the spin current nor the 
charge current are conserved. For a finite repulsive interaction and next nearest neighbor, 
this problem has been studied in the past js?]. One finds two different Luttinger liquids for 
spin up and spin down, which give rise to two different conductances as functions of the spin 
polarization [37] (see equations (26 — 27) and figures (1 — 2) in reference js^]). 
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For the present model [/ = oo, we observe the following: when the Zeeman interaction A 
is small, we do not expect significant changes with respect to the results obtained for A = 0. 

For large values of the Zeeman interaction , ^' ' > 1, the effect of the con- 

straint is negligible. This can be seen from the bias field Ax which induces a space 
dependent oscillation for the constraint (5f^'^'"""(x) in equation (43) . Due to the 
oscillations we obtain, Qf '^'^'"""(x) ~ 0. The wire Hamiltonian is replaced by an 
unconstrained polarized wire H^i^e ~ -f^iire = Ha dx\ ^ ^ ^ -(^x^ffa^))^ + {dx'd^{x))'^] 



Therefore, we recover the robust 0.5 plateau in agreement with 22 |. 
8. Conclusion 

We have solved the problem of exclusion of double occupancy using Dirac's method for 
constraints. We have found that the anomalous commutation rules are modified causing the 
conductance to be anomalous. Applying this theory to quantum wires, we show that our 
theory can explain the anomalous conductance observed by and is in agreement 

with the Monte-Carlo simulation reported in refs. 
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FIG. 1: The Id conducting channel: a) noninter acting leads U = regions > | , strongly 
interacting region U = oo restricted to |x| < | and contact region e. 
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100 " 

FIG. 3: The effective voltage on each side of the wire V^+^ = rae(:^)(^)VK[f , |] where W[ 
is a dimensionless function. 
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